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Abstract. The nature of the tenuous disks around A stars has raised considerable controversy in the literature 
during the past decade. The debate whether or not the disk around /3 Pictoris contains gaseous molecular hydrogen 
is only the most recent example. Since CO is in general a poor tracer for the gas content of these low mass disks, we 
discuss here detailed emission line calcul ation s for alternative tracers like C and C + , based on recent optically thin 
disk models by Kamp & van Zadelhoff (2001). The [Cn] 157.7 pm line was searched toward Vega and /3 Pictoris 
— the most prominent A stars with disks — using ISO LWS data, and a tentative detection is reported toward 
the latter object. From a comparison with emission line observations as well as absorption line studies of both 
stars, the gas-to-dust ratio is constrained to lie between 0.5 and 9 for f3 Pictoris. For Vega the [C n] observations 
indicate an upper limit of 0.2 M® for the disk gas mass. Predicted line intensities of C + and C are presented for 
a range of models and appear promising species to trace the gas content in the disks around A stars with future 
instrumental capabilities (SOFIA, Herschel, APEX and ALMA). Searches for CO emission should focus on the 
J= 3-2 line. 


Key words. Stars: circumstellar matter - planetary systems: protoplanetary disks - Stars: early-type - Stars: 
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1. Introduction 


A large number (20 - 40 %) of nea rby A stars is sur¬ 
rounded by dust disks (Cheng et al. |l992 , Habing et al. 
2001 ), see for e xampl e (3 Pictoris (20 Myr; Barrado y 
Navascues et al. 1999| ) or HR4796A (8 Myr; Stauffer et 
al. 1995 ). However, little is known about the nature of 
these disks: are they dusty debris disks or do they still 
have dust and gas and are possibly in a protoplanetary 
phase? 

The disks around A stars have been extensively 
searched for gas. CO rotational line observations revealed 
too little gas compared to the detected dust. Upper lim¬ 
its suggest that the gas mass in these disks can be up to 
a factor of 1000 lower than deduced from the typical in¬ 


terstellar gas-to-dust ratio of 100 (Yamashi ta et al. 1993; 
Sa voldin i & Galletta 1991; Zuck erman et al. 1995 : Dent et 
al. 1995 ; Liseau & Artymowicz |l998| ; Liseau 1999 ). This 
has been interpreted that the gas dissipates very rapidly 
from such disks on timescales of a few million years. 
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In a previous paper, Kamp & Bertoldi ( |2000D showed 
that in low mass optically thin circumstellar disks CO is 
likely to be photo-dissociated either by the stellar or by the 
interstellar ultraviolet (UV) radiation. If effective shield¬ 
ing prevents CO from being photo-dissociated, it is mainly 
found in the form of CO ice mantles on the dust grains 
in the cool outer disk regions. This indicates that CO is 
actually a poor tracer of the gas and one should look for 
alternative tracers. A similar conclusion has been reached 
for younger, more massive disks by Thi et al. (|2001b). 


In order to calculate the line emission from different 
atoms and molecules in these disks a realistic modeling 
of the gas temperature is needed. With the low densities 
found in these disks, gas and dust are not in collisional 
equilibrium; hence the gas and dust temperatures divert. 
The computation of the ga s tem perature was performed 
by Kamp & van Zadelhoff ( 2001 ), including the relevant 
heating and cooling processes for the gaseous species. 


Recently, Thi et al. (2001) found hints of molecu¬ 
lar hydrogen in the disk around /3 Pictoris using the 
Short Wavelength Spectrometer (SWS) on board of the 
Infrared Space Observatory (ISO). They tentatively de¬ 
tected the S(0) and S(l) lines and deduced from the line 
ratio an excitation temperature of ~100 K assuming LTE 
(?r cr > 10 3 cm” 3 ). The resulting H 2 mass is ~54 with 
an uncertainty up to a factor of 3. On the other hand, 
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Lecavelier des Etangs et al. ( 2001 ) failed to detect H 2 line 
absorption superposed on the broad Ovi emission dou¬ 
blet at 1035 A using the FUSE satellite. Most recently, 
Olofsson et al. (2001) observed extended emission in the 
gaseous NaiD lines using high resolution long-slit spec¬ 
troscopy. The observed velocity pattern is consistent with 
resonance scattering arising in a disk in Keplerian rota¬ 
tion. The emission can be traced from less than 30 AU to 
distances of at least 140 AU. 

This short summary of the debate on (3 Pictoris shows 
definitely the need to find suitable gas tracers, which can 
resolve the question whether the disks around young A 
stars are debris disks without any gas, or protoplanetary 
disks in which gaseous planet formation may still proceed. 

This paper concentrates on CO, C and C + as possi¬ 
ble gas tracers in these disks. Although H 2 is contained 
in our chemical model, we abstain from a calculation of 
the near infrared molecular lines, since it would involve a 
detailed treatment of the ultraviolet pumping by stellar 
and interstellar photons which is beyond the scope of this 
paper. Moreover we do not calculate the Na emission from 
our models. Exploratory one-dimensional calculations in¬ 
dicate that in these disk models the NaiD pumping is not 
optically thin for the stellar radiation. Hence any radiative 
transfer has to take into account not only the stellar radi¬ 
ation field, but also the local radiation field arising from 
the Na atoms within the disk. 

In the following sections we shortly summarize the ba¬ 
sic properties of the optically thin disk models and de¬ 
scribe the method used to derive the level populations 
and the line emission. We present results for a number 
of disk models with varying parameters and discuss them 
in the light of recent observations. In addition ISO LWS 
observations of the [Cil] 157.7 p m line toward Vega and 
/3 Pictoris are presented. At the end, suitable tracers to 
detect the gas in disks around A stars are discussed, and 
line intensities are calculated for future instruments and 
facilities. 


2. The disk models 

We use here the models for low mass disks surround¬ 
ing A stars with stellar spectra similar to (3 Pictoris 
and Vega presented in an earlier paper (Kamp & van 
Zadelhoff 2001). We shortly summarize the main fea¬ 
tures of the models and refer to the two previous papers, 
Kamp & Bert oldi (|2000| ; hereafter Paper i) and Kamp & 
van Zadelhoff ([2001; hereafter Paper n), for further details. 


2.1. Basic physics of the models 

We assume thin hydrostatic equilibrium models 

n(r, z) = tl\ (r/Ri)~ 


-2.5 £ -z /2h 


( 1 ) 


surface density is —1.5, in reasonable agreement with the 


literature values of brightness profiles (Hayashi et al. 1985 
Dutrey et al. 1996|: Augereau et al. [2001 ) ranging from — 1 
to —1.5. 

The radiation field F v is composed of a stellar and 
an interstellar component. The stellar field is described 
by an ATLAS9 photospheric model (Kurucz 1992) for the 
appropriate stellar parameters, while the interstellar com¬ 
ponent is described by a Habing (1968) field with a flux of 
Fj-i = 1.2xl0 7 cm -2 s -1 from 912 A to 1110 A penetrating 
homogeneously through the entire disk. 

The dust temperature follows from radiative equilib¬ 
rium assuming large spherical black body grains of size a 


Tdust = 282.5 (L*/L 0 ) 1/5 (r/AU)” 2/5 {a/pm)~ 1/5 , (2) 

with the stellar luminosity in units of the solar luminosity 
Lq. The assumption of radiative equilibrium is correct 
for the optically thin disks described in this paper. For 
higher mass, optically thick disks, the dust-temperature 
will depend on a detailed calculation of the transfer of 
photons, including scattering, through the disk. Scattering 
will in that case depend strongly on the dust-grain size. 
The use of a grain size distribution would directly affect 
the dust temperature and the shielding of UV radiation; 
the latter is negligible for the tenuous disks discussed in 
this paper. The dust temperature enters the chemistry 
e.g. via the formation of H 2 and freezing out of CO and 
influences the gas temperature in some parts of the disk 
via IR pumping of e.g. OI fine structure lines. Moreover 
the photoelectric heating of the gas depends on the dust 
grain size distribution. For the purposes of this work, the 
grain size distribution is approximated by an effective dust 
grain size, which represents the mean properties of the 
dust phase. Hence, in the following a single grain size of 
3 pm is assumed for all calculations. 

The gas temperature and the chemical composition of 
the circumstellar disk are derived by coupling the heating 
and cooling balance and the equilibrium chemistry. The 
most relevant heating processes are: heating by infrared 
background photons, heating due to the drift velocity of 
dust grains, cosmic ray heating, photo-electric heating, 
and heating due to formation and photo-dissociation of 
H 2 . Cooling, on the other hand, is dominated by [Oi] and 
[C Ii] fine-structure lines as well as the rotational lines of 
H 2 and CO. The most uncertain process is the drift ve¬ 
locity heating: as grains are accelerated by the stellar ra¬ 
diation field and moving through the gas they loose part 
of their momentum. The final drift velocities reached by 
the dust grains depend strongly on the rotation of the 
disk — this cancels a certain fraction of gravity — and 
on the amount of momentum transfer assumed. In order 
to bracket reality we obtained disk models for the two ex¬ 
treme cases: no drift velocity and a maximum drift velocity 


with a dimensionless scale-height Ff = h/r = 0.15. The 
inner and outer radius of the disk are fixed to R\ = 40 AU 
and R 0 = 500 AU. The power-law exponent of the disk 


max _ 
"drift — 


l(U 


2 - 1 - 7) 4 

rad ' u g a 


x 0.5 


0.5 


cm s 


( 3 ) 
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where v gas is the thermal velocity of the gas, and / ra( j is 
the radiation pressure on the grains. 

The chemical network consists of 47 atomic, ionic, and 
molecular species that are related through 260 gas-phase 
chemical and photoreactions. A number of reactions is 
treated in more detail like H 2 and CO photo-dissociation, 
and C ionization. The only surface reactions incorporated 
are H 2 formation and freezing out of CO on cold dust-grain 
surfaces. Since we are dealing with large dust particles, we 
reduced the H 2 formation rate according to the reduced 
grain surface area. The abundance of CO ice is due to a 
balance between freezing out of gaseous CO and reevapo¬ 
ration of CO ice. A modified Newton-Raphson algorithm 
is used to obtain a stationary solution of the entire chem¬ 
ical network. 


2.2. Grid of A star disk models 

We chose in Paper 1 and 11 two very prominent represen¬ 
tatives of the class of A stars with disks around them: 
Vega and /? Pictoris. They differ mainly with respect to 
the strength of their radiation field, with Vega, spectral 
type AOV, having an integrated ultraviolet flux about 
6 000 times larger than the interstellar UV field at 40 AU, 
while the /3 Pictoris UV field, spectral type A5V, is weaker 
than the interstellar held outwards of 53 AU. At 40 AU it 
is about 1.8 times the Habing held. These two radiation 
helds roughly give the lower and upper limits for A stars. 


Table 1. Overview of the radiative transfer calculations. 
For the models indicated by the lines of C, C + , and CO 
were calculated. 


F v 

T?*a p* 

J - 1 / ± u 

Ft 

1 riS 
x V 1 x V 

Ft + F “ 

^drift 

„ max 
^drift 

0 

„.max 

Adrift 

0 

Tgas 

Tdust a = r 6 

A = r 

A = r 

A = r 


Star M e 


[M®] 


A5V 0.2 
A5V 2.0 




AOV 0.2 
AOV 2.0 
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V 


V 

V 


V 

V 


V 

V 


Note: a F* and F^ s denote the stellar and interstellar radiation 
held, respectively 

b A = r denotes that the gas temperature is derived 
from the heating/cooling balance. 


We use the disk models derived in Paper 11 for both 
types of stars and we add two more models for the 
/3 Pictoris case, where we include an interstellar UV ra¬ 
diation held penetrating homogeneously through the disk 
without any shielding. This is a valid assumption for disk 
masses of 2 and 0.2 M® (see Paper 1 ). Due to the high UV 


hux of an AOV star the interstellar radiation held can be 
neglected for Vega. 

The mass of the disk around (3 Pictoris is 54 M® de¬ 
rived from the H 2 S(0) and S(l) emission lines (Thi et 
al. 2001) and 44 M® from dust emission, assuming a con¬ 
stant gas-to-dust mass ratio of 100 (Chini et al. |1991 ). In 
this paper only calculations for lower disk masses are per¬ 
formed, 0.2 and 2.0 M®, which bracket the observations 
for p Pictoris (upper limit for the CO emission lines and 
Cl column densities from absorption lines). Our models 
require that the heating and cooling processes are optically 
thin. This prerequisite is not fulfilled in the M > 2.0 M® 
disk models since the [0 1 ] cooling lines will become opti¬ 
cally thick. 

In order to show the effects of the inclusion of inter¬ 
stellar radiation as well as to give an overview of the tem¬ 
perature structure and chemical abundances in the models 
presented here, two cases are described in more detail. In 
Fig. 0 two disks with a mass of 2.0 M® are shown irra¬ 
diated by a star of spectral type A5V. The disk model 
without interstellar UV radiation field (right hand side of 
Fig. [I]) was presented in Paper II. The new model includes 
now also the interstellar UV field (left hand side of Fig. 
[j]). For each model the temperature structure of the gas 
and dust are given as well as the densities of the molecu¬ 
lar, atomic and ionic species at the midplane and at one 
scale height (h = 0.15r) of the disk. Both models are cal¬ 
culated under the assumption that Udrift = 0. To show the 
effects of the drift velocity on the gas and dust temper¬ 
atures, the corresponding models with Udrift = ^drift are 
presented as well. As the chemistry depends only weakly 
on the temperature, at least in the temperature range cov¬ 
ered in these models, the densities of each of the species 
_are very similar with or without taking Vdrift into account. 

In the model without interstellar UV radiation field, 
hydrogen is mostly molecular and carbon is in the form of 
_CO. As the dust temperature drops below 50 K, the crit¬ 
ical value assumed in these models for CO molecules to 
freeze out onto grain surfaces, CO ice becomes the dom¬ 
inant carbon bearing species. The abundances of C and 
C+ drop according to the freeze out of the gaseous CO. 
The inner disk region has gas temperatures around 10 K 
jiue to efficient CO rotational line cooling. Outwards of 
~ 200 AU, 0 1 fine structure fine cooling takes over, lead¬ 
ing to slightly higher gas temperatures, T gas ~ 20 K. The 
.main heating processes are pumping of the 0 1 fine struc¬ 
ture levels by infrared background photons (inner part of 
the disk) and cosmic ray heating. If we include drift veloc¬ 
ity heating, it becomes the dominant heating mechanism 
and as the cooling processes remain the same, the equilib¬ 
rium gas temperature is higher. Even though the chem¬ 
istry hardly changes due to the higher gas temperature, 
the excitation of atoms and molecules depends strongly on 
T gas . Therefore, detailed modeling of emission lines from 
circumstellar disks has to be based on disk models with a 
realistic determination of the gas temperature. 

The main difference arising from the inclusion of an in¬ 
terstellar UV radiation field is the enhanced dissociation 
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Fig. 1 . Example of two of the modeled disks with 2 and stellar UV radiation appropriate for a /3-Pictoris (A5V) 
star. The panels (a)-(d) show a disk with interstellar UV field, the panels (e)-(h) the same disk without the interstellar 
UV radiation field. From top to bottom the panels represent: (a) and (e) the temperature of the gas (greyscale) and 
the dust (white contour) for a disk with Adrift = 0; (b) and (f) the density of the different species at the midplane; (c) 
and (g) the density of the different species at 1 scale height (h = 0.15r); (d) and (h) temperature of the gas (greyscale) 
and the dust (white contour) when Udrift = ^drift- The adopted lines in (b), (c), (f) and (g) are: H: solid, H 2 : dotted, C: 
dashed, C + : dash-dotted, CO: dash-dot-dotted, COi ce : clash-dot-dotted (thick), e-: long dashed, T<j U st: dashed (thick) 
and T gas : solid (thick). 


of H 2 and CO and the enhanced ionisation of C. The CO 
ice abundance (CO; ce ) is only marginally affected. In the 
disk with interstellar UV, CO is constantly formed and 
dissociated, but most of the CO freezes immediately out 
onto the grains. The CO; ce is the sink for all C-bearing 
species at larger radii, except for the low density surface 
layer, where hardly any CO freezes out. Since the initial 
ratio of [0]/[C]> 1, O together with C + and H 2 are the 
most important coolants available. Due to the enhanced 
UV flux, H 2 dissociation and formation are the most im¬ 
portant heating sources, leading to higher equilibrium gas 
temperatures than in the previous case without the in¬ 
terstellar UV radiation field. As soon as C + is lost as a 


coolant because of its incorporation into CO ice, the gas 
temperatures rise quickly to values in excess of 100 K. 
Incorporating the drift velocity into this model does not 
alter the gas temperature, because H 2 dissociation and 
formation heating are more efficient than drift velocity 
heating. 

A complete list of all adopted models and performed 
radiative transfer calculations can be found in Table [I]. 


3. Radiative transfer calculations 


For the line radiative transfer the Monte Carlo code de¬ 
veloped by Hogerheijde & van der Tak ( [2000 ) is used. 
The disk models described in Sect. 0 are interpolated on 
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a cylindrical grid (26 x 10 grid cells), with a logarith¬ 
mic radial grid in order to smoothly follow the density 
and temperature gradients. The size of the grid is small 
due to computational constraints. A comparison run with 
twice the number of grid cells gave a less than 20% differ¬ 
ence in integrated line-intensity, comparable to the obser¬ 
vational errors. The calculations concentrate on the non 
local-thermal-equilibrium (NLTE) intensities of the fine 
structure lines of C and C + and the rotational transitions 
of CO. Table shows the respective line data and beam 
sizes used for convolution in the section on comparison to 
observations (Sec. j(|). The following paragraphs explain 
the calculations for the individual species in more detail. 

Table 2. Adopted line data for the fine-structure lines of 
C and C+, and the rotational lines of CO. 



Line 


Aui 

[s- 1 ] 

Beam 

["] 

c 

809.3 GHz 

3 P 2 - 3 Pi 

2.65 10” 7 

6“;8 6 ; 

26 c 


492.2 GHz 

3 Pi- 3 Po 

7.93 10” 8 

10“;13 i ’; 

44“ 

c+ 

157.74 pm 

“P3/2- 2 P1/2 

2.29 10~ 6 

ll c ;16 d 

CO 

115.3 GHz 

1-0 

7.17 10" 8 

43“ 


230.5 GHz 

2-1 

6.87 10~ 7 

22“ 


345.8 GHz 

3-2 

2.48 10“ 6 

14“ 


“15 m telescope (SEST, JCMT), 6 12 m telescope (APEX), 
C HIFI (Herschel), d SOFIA 


In the above mentioned code, the equations of statis¬ 
tical equilibrium are solved in an iterative fashion, where 
all photons start at the outer boundary with an intensity 
given by the 2.728 K cosmic background radiation. The 
abundances of the trace species under consideration as 
well as their most important collision partners are taken 
from our stationary disk models (Paper II and Fig. [[]). 

The code allows the use of two different collision part¬ 
ners and Table || shows the most important partners for 
each species and the respective references for the colli- 
sional rate coefficients. Stimulated absorption and emis¬ 
sion by far-infrared radiation from dust is not significant 
for these species. For C-H 2 and C + -H 2 collisions, we as¬ 
sume that the ortho-to-para ratio of molecular hydrogen 
is determined by the local gas temperature. In the case of 
C+, scaled H 2 collision cross sections for atomic hydrogen 
are used; hence the total density of the collision partner 
is given by the weighted sum of the H 2 and H density, 
namely n(H 2 ) +n(H)/0.57, and n(e) is considered as the 
second partner. Comparing CO-H and CO-e _ collision 
cross sections, it is found that for n(H)/n(e“) < 10 3 , col¬ 
lisions with electrons are more important than collisions 
with neutral hydrogen. In the case of all AOV models and 
the A5V models with interstellar radiation, n(H) is the 
main collision partner for CO; in the remaining models 


Table 3. Collision partners for the different species and 
the respective references for the collisional rate coeffi¬ 
cients. 


Species 

Collision partner 

Reference 




C 

H 2 

Schroder et al. 

(1993) 

C + 

H 

H 2 

H 

Launay 
Flower ( 
Flower ( 

fc R 

1977 

1977 

aueff (1977) 

) 


e“ 

Mendoza (198E 

L 


CO 

H 2 

Keenan et al. ( 
Schinke et al. ( 

1986 

1985 

) 

) 


H 

e - 

Chu & Dalgarno (1975) 


“Scaled from H 2 collision rates. See text for details. 


not enough H 2 is dissociated and collisions with electrons 
dominate over atomic hydrogen collisions. 

The resulting populations at each position in the disk 
are used to compute the sub-millimeter line profiles of 
CO, C, and C + using a ray tracing program which cal¬ 
culates the sky brightness distribution. Observations of 
circumstellar ultraviolet absorption lines of CO and C 
around (3 Pictoris reveal a li ne-br oadening parameter b = 
1.3 km s _1 (Roberge et al. |2000| ). This value is used for 
all model calculations. 

4. Results 

As a first result the line profiles are calculated for a beam 
with the same size as the apparent diameter of the disk on 
the sky. The beam is centered on the stars and has a size 
of 50.5" for simulating (3 Pictoris and 128.9" for Vega at 
their appropriate distances, 19.82 and 7.76 parsec, respec¬ 
tively. The adopted inclinations are taken from (3 Pictoris 
(edge on) and Vega (face on). 

Since the emission lines are optically thin, the inte¬ 
grated intensities can be converted directly to any incli¬ 
nation. These results can therefore be used for any disk 
around an A star as long as the total integrated flux 
f T m bdv x is kept constant. In the following the pre¬ 
sentation of the results is separated into two sections. The 
first deals with A5V and the second with AOV stars. 

4.1. Spectral type A5V 

We choose atomic carbon as an example to illustrate the 
effects of the various model parameters on the resulting 
line profiles (Fig. |2j). The shapes of the emission-line pro¬ 
files show a double peaked structure, with subtle differ¬ 
ences due to excitation and/or abundances (e.g., Figs. |2f 
and |^b) . 

For the 2M ffi model including the interstellar UV radi¬ 
ation field (Fig. ||c,d), the gas temperature rises above 
50 K resulting in a 3 P 2 - 3 Pi/ 3 Pi~ 3 Po ratio of approx¬ 
imately unity. The carbon density rises by an order of 
magnitude compared to the other 2M® models, because 
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Fig. 2. [Cl] fine-structure line profiles (main beam bright¬ 
ness temperature versus rest velocity) in different A5V 
models: (a) 0.2 M® model with Wdrift = and an in¬ 
terstellar radiation field, (b) 0.2 M® model with thrift = 0 
and interstellar radiation field, (c) 2 M® model with 
Adrift = ^drift and interstellar radiation field, (d) 2 M® 
model with Udrift = 0 and interstellar radiation field, (e) 
2 M® model with Udrift = ^drifti (f) 2 M® model with 
V drift = o, (g) 2 M® model with T gas = T dust . The solid line 
denotes the 3 Pf 3 Po transition at 492.2 GHz (609.13 /jm) 
and the dashed line the 3 P 2 - 3 Pi transition at 809.3 GHz 
(370.42 /im.) 


CO is no longer the dominant reservoir of carbon due to 
the additional photodissociation. There is no difference in 
gas temperature between the calculation with or without 
drift velocity heating, as the temperature balance is in this 
case entirely determined by Ho formation and dissociation 

(Fig-0)- 

For the cases with only the stellar radiation field, the 
difference between the 2 M® models with (Fig. |°) and 
without (Fig. gf) drift velocity is due to the gas tempera¬ 
ture. T gas drops below 20 K (Fig. [I]) for a large part of the 
disk without any drift velocity heating, keeping most of 
the carbon in the ground state. Hence the lines are signif¬ 
icantly weaker than with drift velocity heating included. 

In the case where T gas = T dust (Fig. gg), the 3 P 2 - 
3 Pi line is stronger than the 3 Pi- 3 Po line contrary to the 
other 2 M® models (Figs, go, f) which include heating and 


cooling. This is a direct result of the temperature. In the 
models with heating and cooling, the gas temperature is 
well below 30 K in the inner 150 AU, while it stays above 
50 K if T gas = T dust is assumed. The 3 P 2 fine-structure 
level of C lies at 62 K and is only significantly populated 
at temperatures above 30 K. 




-10 -5 0 5 10 -10 -5 0 5 10 

v [km/s] v [km/s] 


Fig. 3. [C n] 157.74 fj ,m fine-structure line profiles (main 
beam brightness temperature versus rest velocity) in dif¬ 
ferent A5V star models: (a) 2 M® model with Vdrift = '^drift 
and interstellar radiation field, (b) 0.2 M® model with 
Adrift = ^dSt and interstellar radiation field, (c) 2 M® 
model with Udrift = and (d) 2 M® model with 

Adrift = 0 (note the flux is multiplied by 10 3 ). 


If the disk mass is reduced by a factor of 10 (Figs, gi, 
b), the carbon density decreases by a slightly larger fac¬ 
tor, because the ionization balance shifts towards C + as 
shielding decreases in lower mass models. The tempera¬ 
ture stays above 40 K in the inner 100 AU of the disk, 
giving rise to equally strong emission lines. This is not the 
case if the drift-velocity heating is absent. In that case the 
temperature drops below 30 K causing the 3 P 2 - 3 Pi line 
to become much weaker and the 3 Pr 3 Po slightly higher. 

For C+ and CO we concentrate on four models, the 
2 and 0.2 M® model with interstellar UV radiation and 
Adrift = '(’“rift; and th e 2 M® models with u drift = thrift and 
Adrift = 0 without interstellar radiation. The combination 
of these models shows all the possible influences in the 
models. Comparisons between panels (a) and (c) show the 
influence of the interstellar UV radiation field, between (a) 
and (b) of the disk mass, and between (c) and (d) of the 
drift-velocity heating on the emission lines of C + and CO. 

The influence of the interstellar radiation field on the 
temperature and chemical abundances can be seen in 
Fig. [T[ The CO abundance is reduced due to enhanced 
photodissociation, and the C + line therefore becomes 
much stronger. Moreover, the lower CO transitions J=l- 
0 and J=2-l are generally weaker than the J=3-2 line. 
This is due to typical gas temperatures in excess of 30 K 
(the <7=3 level lies at 33 K) due to the interstellar ra¬ 
diation field. The level populations were calculated using 
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Fig. 4. CO rotational line profiles (main beam brightness 
temperature versus rest velocity) in different A5V star 
models: (a) 2 M 0 model with Udrift = an d interstellar 
radiation field, (b) 0.2 M 0 model with thrift = and 
interstellar radiation field (note the flux is multiplied by 
10 2 ), (c) 2 M 0 model with v dr ift = and (d) 2 M 0 

model with u dr ift = 0. Solid line: J=l-0, dashed line: J=2- 
1, dot-dashed line: J=3-2. 


Fig. 5. [Cl] fine-structure line profiles (main beam bright¬ 
ness temperature versus rest velocity) in different A0V 
models: (a) 0.2 M 0 model with Udrift = ^drift> 0 3 ) 0.2 M 0 
model with u dr ift = 0, (c) 2 M 0 model with u dr ift = 

(d) 2 M 0 model with Udrift = 0. The solid line denotes the 
3pi 3p o transition at 609.13 pm and the dashed line the 
3 P 2 ~ 3 Pi transition at 370.42 /im. Note that the fluxes in 
panels (a) and (b) are multiplied by 10 2 . 


the full statistical equilibrium calculation. This has only 
a marginal effect on the three lowest levels of CO for the 
2.0 M 0 model. Only in the surface layers of the disk are 
the densities low enough to prevent collisionally excita¬ 
tion of J=3. These slight differences are therefore only 
marginally visible in the line-profiles and not distinguish¬ 
able with current telescope facilities. 

The influence of the drift-velocity on the line emission 
is large in the case without interstellar UV radiation. Due 
to the extremely low gas-temperatures in the inner regions 
of the disk, where the CO density is highest, the J=1 level 
is excited more efficiently than the higher levels, making 
the lowest transition the strongest (Fig. lO- 

In the case of the lower mass (0.2 M 0 ) disk, the level 
populations do show clear NLTE effects on the level pop¬ 
ulations, even close to the midplane. The total emission, 
however, is extremely small. Whereas the C and C + line 
emission scale roughly with mass, the CO line emission is 
reduced by nearly three orders of magnitude due to en¬ 
hanced dissociation in the lower mass disk. 

4.2. Spectral type A0V 

Since we adopted the disks around A0V stars to be seen 
face on, no coherent velocities are seen in the line of sight 
and only the micro-turbulent (b = 0.91 km s -1 ) and ther¬ 
mal velocities remain. 

All line profiles show a single gaussian emission peak, 
and are generally orders of magnitude weaker for C and 
CO than in the A5V models. Due to the high UV flux 
from the A0V stars, most of the CO in the gas phase 
is dissociated. The atomic carbon is subsequently ionized 


enhancing the C + abundance and its emission by a few 
orders of magnitude. 

Similar to the A5V stars, we illustrate the differences 
arising from the various models by using C as an exam¬ 
ple. The models with drift-velocity heating are generally 
warmer than those without. This is nicely illustrated in 
Fig. where the 3 P 2 - 3 Pi/ 3 Pi- 3 Po line ratio is higher 
for the models with Udrift = ^drift- The temperature effect 
is more pronounced in the 0.2 M 0 model, as the ratio is 
larger than unity. For the Udrift = 0 models, the cooling is 
more efficient in the low-mass disk giving a relatively large 
3 Pi 3 P 0 line emission compared to the higher transition. 
In the low mass models, the ionization fraction is much 
higher due to a general lack of shielding. This leads to C 
emission that is a factor of 100 lower than in the 2 M 0 
models. 

For C+ the differences are entirely due to the excita¬ 
tion. The C+ densities scale with the mass, but the popu¬ 
lation of the upper level (91.2 K) depends strongly on the 
gas temperature. For CO, the J= 3-2 line is the strongest 
rotational transition in the models with Udrift = ^drift aud 
roughly equal to the 2-1 line in the other two models. 

4.3. Integrated intensities 

The line profiles presented in the above sections are cal¬ 
culated for a single position using different beam-sizes. 
Therefore they do not reveal in which part of the disk the 
emission arises, nor at which position the beam should be 
centered to pick up the maximum emission from the disk. 
To illustrate the radial dependence of the emission, the 
disk is scanned with a small beam over the disk midplane, 
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Fig. 6. [Cii] 157.74 /mr fine-structure line profiles (main 
beam brightness temperature versus rest velocity) in dif¬ 
ferent AOV star models: (a) 0.2 M® model with thrift = 
(b) 0.2 M® model with i>drift = 0, (c) 2 M® model 
with Adrift = Adrift (note that the flux is multiplied by 0.1), 
and (d) 2 M® model with Adrift = 0. 





Fig. 7. CO rotational line profiles (main beam brightness 
temperature versus rest velocity) in different AOV star 
models: (a) 0.2 M® model with Udrift = VdJfft, (b) 0.2 M® 
model with Vdrift = 0, (c) 2 M® model with thrift = ^drift; 
and (d) 2 M® model with Udrift = 0. The solid line is J=l- 
0, dashed line J=2-l, dot-dashed line J=3-2. Note that 
the fluxes are multiplied by the factors indicated. 


2 " in the case of the edge-on A5 V star and 5 " in the case 
of the pole-on AOV star. 

Fig. shows that the maximum emission for the A5V 
star for all three tracers, CO, C and C + is at 2.5", corre¬ 
sponding to 50 AU at a distance of 19.8 pc. In the inner 
disk the particle densities are highest whereas outwards 
the density drops with r~ 2 5 . Beyond 10", i.e. 200 AU, 
carbon is mainly in the form of CO ice (see Fig. [j]) and 
hence the emission gradually drops off. 

In the case of the A0 V star the stellar UV flux is strong 
enough to keep C + as the main emitting carbon-species 
throughout the disk. The emission is strongest between 5 
and 10" (40-80 AU) for C+ and C and between 10-15" 



Fig. 8. Integrated flux along the disk midplane of the 
2 M® model illuminated by an A5 V star, with thrift = 0 
and interstellar radiation field using a 2" beam size: CO 
3-2 (solid line), [Cl] 3 Pi- 3 Po (dashed line), and [CII] 
(157.74 fin i, dash-dotted line). The distance to the star 
is taken to be 19.8 pc. 


(80-120) for CO. Closer to the star the beam is diluted 
due to the inner hole of 40 AU, hence the rapid decline in 
emission close to the star. 



r [arcsec] 


Fig. 9. Integrated flux along the disk radius of the AOV 
star, 2 M® model with Adrift = 0 with a 5" beam size: 
CO 3-2 (solid line), [Cl] 3 Pi- 3 Po (dashed line), and [Cii] 
(157.74 fin i dash-dotted line). The distance to the star is 
taken to be 7.8 pc. 
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5. [Cn] observations of Vega and /3 Pictoris 


The Infrared Space Observatoty (ISO) data archive was 
used to retrieve all observations made with the Long 
Wavelength Spetrometer (LWS) in a radius of 15' to¬ 
ward Vega and f3 Pictoris. This yielded data from projects 
by MBARLOW and RSTARK taken in full grating scan 
mode (AOT L01) and grating range scan mode (LWS 
AOT L02). The spectral sampling interval of all data from 
these projects is 1/4 of a resolution element. The number 
of grating scans (forward and backward) varied between 6 
and 24. Each scan yields simultaneous data accross ten de¬ 
tectors. The total integration time varied between 1172 s 
and 3592 s per AOT. We only considered the data from 
detector LW04 whose central wavelength is at ~ 160/im 
(grating in nominal position). The effective beam size of 
this detector is 78" (FWHM), the spectral resolution is 
about 0.6 /im which corresponds to a velocity resolution 
of about 1150 km s -1 at 157.7 /jm (ISO Handbook Vol. 
IV; Version 1.2). 


The data were automatically processed through the 
LWS pipeline (OLP version 10), yielding LWS Auto anal¬ 
ysis (LSAN) files which contain the flux and wavelength 
calibrated spectrum of an AOT. The absolute f lux ca libra- 
tion is estimated at about 30% (Swinyard et al. 1996 ). The 
wavelength accuracy for the LW detectors is ~ 0.15 /im 
(ISO Handbook Vol. IV; Version 1.2). Final data reduc¬ 
tion was done manually using the ISO spectral Analysis 
Package (ISAP) Version 2.1. Glitches and their related 
decays due to Cosmic Rays were removed from each indi¬ 
vidual grating scan. A medium clip was subsequenly ap¬ 
plied to all scans for all data in a bin which is more than 
2.5 a larger or smaller than the median. Since the observed 
sources are not bright nor larger than the LWS beam, no 
defringing was applied. Subsequently an average to the 
mean of all subscans was done for each bin. A linear base¬ 
line was fitted to the averaged spectrum and the line fluxes 
were estimated through fitting of a gaussian profile with 
a FWHM comparable to the instrumental resolution. 


The resulting spectra are presented in Fig. [h| The 
spectrum toward Vega does not reveal any spectral line, 
a 2 a upper limit to the [C n] line intensity was estimated 
at < 1.6 x 10 -20 W cm -2 pm -1 in a 0.15 /inr bin. The 
line-of-sight toward (3 Pictoris shows a 4 cr feature with a 
maximum intensity of 1.8 x 10~ 20 W cm -2 /inr -1 and an 
integrated intensity of about 1 x 10~ 20 W cm -2 centred at 
157.85 /im, which agrees with the wavelength of the [Cn] 
line within the LWS accuracy. The off-spectrum does not 
reveal any emission at this frequency with a 2 a upper 
limit of 1.0 x 10 -20 W cm -2 ^m _1 in a 0.25 (im bin which 
indicates that there is no contamination by an extended 
interstellar [Cii] emission component. However, a contri¬ 
bution of weak interstellar [C II] cannot be fully excluded 
since the integration time used for the OFF position is a 
factor 2.5 lower than used toward [3 Pictoris. 



155 160 

X(gm) 


Fig. 10. ISO LWS spectra of the [Cn] 157.7 /im line 
for (3 Pictoris (ON 05h47ml7.11s -51d03’59.5” and OFF 
051r47m22.66s -50d57’28.5”, J2000). 


6. Comparison of models with observations 


The attempts to detect the gas component in low mass 
disks around A stars concentrate mostly on prominent A 
stars like (3 Pictoris and Vega. In order to compare the 
models to the observations, we have chosen the same beam 
sizes as used in the different observations. Table shows 
the beam sizes for each line and telescope, while Fig. [ll| 
illustrates the areas of the disk contained in the different 
beams for /3 Pictoris (edge-on at a distance of 19.8 pc) and 
Vega (pole-on at a distance of 7.8 pc). For Vega a radial 
offset of 14.7" (114 AU) is applied because the smallest 
beams contain only the star and the inner hole o f the 
disk model. Observations by Holland et a l. ( 1998 ) and 
subsequent modeling by Dent et al. (200C) substantiate 
the presence of these holes. 


6.1. (3-Pictoris 

For (3 Pictoris, the CO J=l-0 and 2-1 lines have 
been searched several times with the SEST (Savoldini & 
Galletta 1994; Liseau & Artymowicz 1998) and the JCMT 
(Dent et al. 1995| ). In all cases only an upper limit on 
the CO column density was derived, the most strict being 
Nco < 3 10 14 cm -2 so far. This is deduced from the CO 


2-1 upper limit found by Liseau & Artymowicz (1998), 
where they obtain a noise level of T rms = 11 ink using 
a binning of 0.9 km s _1 . With a binning of 5 km s _1 , 
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Fig. 11. Density distributions for (3 Pictoris (top) and 
Vega (bottom). The telescope beams used in the calcu¬ 
lations are superposed as dashed circles. See Table || for 
the corresponding telescopes. The beam sizes in the lower 
panel correspond to those in the top panel. 


which corresponds to the halfwidth of the CO line, the 
noise level goes down to 4 mK. To compare these limits 
to the models, T rms has to be converted using the main 
beam efficiency ? 7 m b = 0.38 of the SEST at the corre¬ 
sponding wavelength. The models presented show that the 
,7=2-1 transition in the (3 Pictoris model with 2 and 
the interstellar radiation field reaches a flux of 12.2 mK 
with a binning of 5 km s -1 , while the <7=1-0 line reaches 
only 0.9 mK (Fig. Ha)- At the typical disk temperatures 
of 50 K in the models, mostly the J=2 and <7=3 levels 
are populated. The predicted <7=2-1 line would not have 
been detected by the latest SEST observations, which have 
T m b ,rms — Trms/^lmb — 10.5 UlK. 

The models show that in fact the most favorable line 
to observe is the ,7=3-2 line. A spectrum for this line, 
taken from the JCMT public archive, has a noise level of 
T r ms = 27 mK after smoothing to 1.1 km s -1 velocity 
bins. A much deeper integration should be possible with 
the dual polarization B3 receiver. The current limit is a 


factor of two lower than the modeled line (Fig. [l^a, dot- 
dashed line). 
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Fig. 12. CO rotational line profiles (main beam brightness 
temperature versus rest velocity) in different (3 Pictoris 
models convolved with the observational beams: (a) 2 M® 
model with Adrift = ^drift an d interstellar radiation field, 
(b) 0.2 M® model with thrift = ^dna an d interstellar ra¬ 
diation field, (c) 2 M® model with Vdrift = with¬ 

out interstellar radiation field, and (d) 2 M® model with 
Adrift = 0. The solid line is ,7=1-0 (43"), the dashed line 
J=2-1 (21"), and the dot-dashed line J— 3-2 (14"). Note 
that the fluxes are multiplied by the factors indicated. 


The 2 M® disk models have significantly lower mass 
than the value of ~ 50M® inferred from the dust using 
a gas/dust ratio of 100:1. Scaling the mass up by that 
factor (x25), both the ,7=2-1 and J=3-2 line should have 
been easily detected even though the CO gas abundance 
is unlikely to scale linearly due to enhanced freeze-out. 

Comparison between Figs. |l2| c and d reveals that with¬ 
out drift velocity heating the disk becomes much cooler 
and the lines weaker. Even the line ratios are affected, in 
the sense that the J=2 1 and <7=3-2 lines are now equally 
strong. Moreover, the 0.2 M® model gives CO lines below 
0.1 mK, indicating that the threshold for CO detection 
occurs between 0.2 and 4 M® (see Fig H>, assuming that 
CO scales linearly in a narrow range of disk mass) for the 
distance and inclination of /3-Pictoris. Roughly a factor 
of 10 in flux is generally gained compared to the aver¬ 
aged intensity from the entire disk (Fig. ||). A comparison 
of Fig. [|d and Fig. |T^ d reveals that the smaller beams 
pick up mostly the warm CO in the inner disk regions. 
In addition the difference in beam dilution is largest for 
the J=2 1 and <7=3-2 lines compared to the entire disk 
beam. 

On the other hand, there are UV absorption-line stud¬ 
ies along the line of sight towards (3 Pictoris. Roberge et 
al. ( [2000 ) observed CI and CO using the HST STIS high- 
resolution echelle spectrograph. Contrary to previously de¬ 
tected Cl absorption (Jolly et al. 199S), an unsaturated 
spin-forbidden transition was observed. This allowed for 
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an improved determination of the ground state C( 3 P) col¬ 
umn of (2-4) 10 16 cm -2 . The CO absorption results indi¬ 
cate a column of (6.3 ± 0.3) 10 14 cm -2 . 

In the 2 M® disk models, the CO column density 
is 810 15 cm -2 , roughly 10 times higher than the ob¬ 
served value. The CI column density from the model, 
~ 110 1 ' cm” 2 , is also a factor 2.5-5 higher than observed 
in the absorption studies. In case of the 0.2 M® disk mod¬ 
els a CO and CI column density of 8 10 12 and 5 10 15 cm -2 
is reached respectively. Both are lower than the observed 
values, leading to the conclusion that the actual mass lies 
in between these models. 


The detected [CII] line provides a reliable estimate of 
the amount of ionized gas in the disk around /3-Pictoris. 
The emission line was modeled adopting the ISO beam of 
78". The integrated emission for the 2 and 0.2 M® disk 
models is 1.9 x 10” 21 and 2.6 x 10” 22 W cm -2 respec¬ 
tively. Both are lower than the observed value of 1 x 10” 2 ° 
W cm” 2 . At first sight, this would indicate that according 
to the model more gas is present in the disk around (3- 
Pictoris. However, this conclusion depends heavily on the 
assumed UV radiation and further modeling is needed. 


Freudling et al. (1995) derived from the non-detection 
of the HI 21 cm line an upper limit of 2 — 5 x 10 19 cm -2 
for the neutral hydrogen column density. The column cal¬ 
culated for the 2 M® disk model including the interstellar 
field is 8 x 10 2 ° cm -2 and for the 0.2 M® disk model 
8 x 10 19 cm” 2 , a few times higher than observed. 


The Hi column density in the 2M® model is cl ose to 
IV(H) = 5x 10 2 ° cm” 2 deduced by Olofsson et al. ( [2001 ) 
from their sodium observations using a solar Na/H abun¬ 
dance ratio. In the models described here, the available 
HI in the gas is directly related to the amount of H 2 self¬ 
shielding of the stellar radiation. The interstellar radiation 
field penetrating through the disk is not shielded by H 2 , 
hence comparable columns of H 2 and H 1 are reached (see 
paper I). A disk with a 25 times higher mass would be 
able to shield the midplane of the disk and the H 1 column 
will decrease. This process would hardly affect the carbon 
species in the outer regions, because CO will freeze out, in¬ 
corporating all the available carbon. In the inner regions, 
the C 1 and CO densities would increase in a higher mass 
disk and the column densities would be too high compared 
with observations. 


Our models show that a /3-Pictoris disk with 2 M® of 
gas might have escaped detection in CO. The absorption 
studies towards [3 Pictoris show that the models give CO 
and C 1 column densities that are up to a factor of 10 too 
high. In case of the 0.2 M® models all of these upper limits 
are met as C + and CO; ce become the main reservoirs of 
carbon. In the case of the H 1 absorption even the low mass 
0.2 M® disk overestimates the observed upper limit. Since 
the estimate by Olofsson et al. (2001) differs by a factor of 
10 and the true H 2 formation rate used is not known, not 
too much weight should be given to the Hi constraint. 
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Fig. 13. CO rotational line profiles (main beam brightness 
temperature versus rest velocity) in different Vega mod¬ 
els, at 14.7" offset and convolved with the observational 
beams: (a) 0.2 M® model with Udrift = 0 3 ) 0.2 M® 

model with thrift = 0, (c) 2 M® model with Adrift = ^drifti 
and (d) 2 M® model with Udrift = 0. The solid line is J= 1- 
0 (43 "), the dashed line J=2-l (21"), and the dot-dashed 
line J= 3-2 (14"). Note that the fluxes are muliplied by 
the factors indicated. 


6.2. Vega 


For Vega, there are in general less observations than for 
/3 Pictoris. Due to the pole-on geometry, absorptio n line 
studies for Vega are impossible. Yamasliita et al. ( 1993 ) 
searched for the CO J=l-0 line with the 45 m Nobayama 
Radio Observatory (NRO). They detected no CO emission 
down to a T rms of 33 mK in 1 km s” 1 bins. Even with the 
bette r sen sitivity of T rms ~ 10 mK of the JCMT, Dent et 
al. ( 1995 ) did not detect any CO J=2-l around Vega. 
The above calculations show that a 2 M® model for Vega 
reaches at most 1 mK for the J= 3-2 transition (Fig. [TTijc). 
The lowest J=l-0 transition has fluxes below 0.1 mK. 
This is entirely due to photodissociation. Simple scaling 
suggests that the lines will still not be detectable with 
current instruments if the disk mass is enlarged by a factor 
of 20. 

The difference seen in comparin g th e models with and 
without drift-velocity heating (Fig. |l3p-d) is entirely due 
to differences in gas temperature. Drift-velocity heating 
is the most important heating source of the gas and, if 
omitted, H 2 formation/dissociation as well as photoelec¬ 
tric heating of the /im-size dust particles remain, leading 
to lower gas temperatures. The CO fluxes are a factor of 
10 higher than with the beam of 128.9" corresponding to 
the entire disk (compare Fig. 0 and |l3|). This is due to 
the low CO densities in the outer regions picked up in 
the large beam. The sometimes different line ratios can be 
explained by the restricted area picked up by the smaller 
telescope beams. 

The models for the Vega disk predict [C 11 ] integrated 
emission lines of 1.910” 20 and 7.910” 22 W cm” 2 for the 2 
and 0.2 M® disks. The former integrated intensity corre- 
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sponds to a peak intensity of 3.4 x 1CV 20 W cm -2 /im _1 
and the latter to 1.4 x 10 -21 W cm -2 p,m -1 . The ISO 
satellite observed Vega in the wavelength range of the [C n] 
157.74 fj,m line (Sect. ||) and obtained an upper limit of 
1.6 x 1(V 21 W cm -2 /im -1 . This would suggest a maxi¬ 
mum amount of gas of ~ 0.2 in the disk around Vega. 


7. Future observations 

We present here [C i] predictions from our models appro¬ 
priate for the new MPIfR/SRON 800 GHz and 460 GHz 
single channel receivers to be installed at the Atacama 
Pathfinder Experiment (APEX) in 2003. APEX is located 
at 5000 nr at Chajnantor in Chili and is perfectly placed 
to observe /3-Pictoris. For Vega, the [Cl] emission line is 
calculated using the JCMT telescope instead of APEX. 
Additionally [C i] and [C ii] lines for the future instrumen¬ 
tation on board of SOFIA and Herschel are calculated. 
The velocity resolution of all the telescopes is set at 1 km 
s -1 to detect the line profiles except for the [C i] 809.3 GHz 
(370 ^m) line, where 5 km s -1 is used. 
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Fig. 14. [C i] fine structure line profiles (main beam 
brightness temperature versus rest velocity): (a) 2 M® 
model for /3 Pictoris with Udrift = 0 and interstellar ra¬ 
diation field appropriate for the APEX (note that the flux 
is muliplied by 0.1), (b) same as (a) but appopriate for 
Herschel, (c) 2 M® model for Vega with thrift = 0 ap¬ 
propriate for JCMT, (d) same as (c) but appropriate for 
Herschel. The solid line denotes the 3 Pi~ 3 Po transition at 
609.13 /im (10, 13, and 44" beam for JCMT, APEX, and 
Herschel respectively) and the dashed line the 3 P 2 3 Pi 
transition at 370.42 /im (6, 8, and 26" beam for JCMT, 
APEX and Herschel, respectively) 


Since detailed sensitivity estimates are not yet avail¬ 
able for APEX, the limits in the appropriate bands are 
calculated for the JCMT. In general, APEX, due to its 
excellent site and stable atmosphere has better sensitivi¬ 
ties. 

The [Cl] 370.42 (im intensity of the 2 M® disk model 
with interstellar UV field for [3 Pictoris yields a flux of 


300 rnK, whereas the 2 M ffi Vega model yields only 15 mK. 
Since the source is at a low zenith angle, the current 
370 fim receiver at the JCMT would require 200 min¬ 
utes of observing time for a 5 km s -1 velocity resolution 
to reach a 3er detection for a 300 mK peak intensity of 
the line. This is clearly not feasible. The expected up¬ 
graded receiver at APEX would reach this limit in about 
20 minutes for /3-Pictoris. A higher spectral resolution of 
1 km s _1 would increase the time to about 8 hours. 

The lowest [Cl] line at 492.2 GHz (609.13 /iin) is cal¬ 
culated using the JCMT specifications. A typical noise 
level of T\= 40 mK in 1 km s _1 bin could be reached 
in 1.1 hours of integration time at an elevation of 60°, 
appropriate for the /3-Pictoris disk from the APEX site. 
The predicted T m b = 100 mK line should be detectable 
with this configuration. For a star like Vega with a strong 
UV radiation field, this line is not observable in low mass 
(< 2 M®) disks. 

For the HIFI instrument on board of the Herschel satel¬ 
lite, the [C i] line at 609.13 //m is at the lower limit of Band 
1, while the 370.42 /jm line falls in Band 3. The sensitiv¬ 
ity estimates for a long integration (5 hours) are 2 and 
3 mK respectively (de Graauw & Helmich 2001). Since 
the Herschel beam is much larger than the JCMT beam, 
nearby edge-on disks will suffer from large beam dilution. 
Still, both the edge-on /3-Pictoris disk and the pole-on disk 
around Vega should be just observable in both lines. 
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Fig. 15. [Cll] 157.74 /rrn fine structure line profiles (main 
beam brightness temperature versus rest velocity): (a) 
2 M® model for (3 Pictoris with Udrift = 0 and interstellar 
radiation field, (b) 2 M® model for Vega with Vdrift = 0. 
The thick line denotes the Herschel beam of 11", the thin 
line the SOFIA beam of 16". Note that the fluxes are 
muliplied by the factors indicated. 


The [Cn] fine structure line at 157.74 (im can be ob¬ 
served in Band 6 of the HIFI instrument on board of 
Herschel. In a long exposure (5 hours), a sensitivity of 
~6 mK is reached (de Graauw & Helmich 2001 ). A 2 M® 
model will be observable in both Vega and /3-Pictoris. 

If the slightly larger beam of SOFIA is used, 16" com¬ 
pared to the 11" of Herschel, the model fluxes for the [C n] 
line decrease by about a factor 2 for an edge-on disk model 
due to additional beam dilution. Since we applied an offset 
of 15 " to the pole-on disk model, there is hardly any dif¬ 
ference between the two beam sizes for this model. In both 
cases, the beam is entirely filled by the disk model. Apart 
from the difference in beam size, HIFI and the GREAT in- 
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strument on SOFIA will have comparable sensitivity lim¬ 
its and integration times. 

The Atacama Large Millimeter Array (ALMA) will 
reach a sensitivity and angular resolution which will be 
orders of magnitude larger than the current status. In the 
case of Vega and (3 Pictoris, angular resolution is of minor 
importance due to the relatively small stellar distances. 
Adopting the expected sensitivities, an hour of integration 
on the CO 3-2 line with a spectral resolution of 0.5 km s -1 
would give a 5 cr detection of 0.05 K. This result is obtained 
for the compact configuration with baselines up to 150 m, 
which gives a spatial resolution of 2 ". Figure |8| shows that 
the integrated line intensity of the CO 3-2 line peaks at 
a 2" offset for /f-Pictoris. Using a spectral resolution of 
0.5 km s -1 the local emission at 2" offset reaches 2.1 K. 
At the wavelength of the [C i] 3 Pi — 3 P 0 transition the 5 cr 
level reached with ALMA is expected to be twice as high 
(0.1 K). The modeled intensity for this line peaks at 1.2 K. 
In both cases the CO and [Cl] line should be detectable 
in a few seconds and half a minute respectively. For Vega 
the detection limits are more stringent and peak intensi¬ 
ties (with a 2" beam) will be 0.05 for the [Cl] line and 
6 x 10 -4 K for the CO line. [Cl] will thus be detectable 
in 4 hours integration time and CO is undetectable (7000 
hrs) for this binsize at a 5 a level. 


8. Discussion and conclusion 


The optically thin models described in this paper provide 
a tool to constrain the gas mass in circumstellar disks on 
the basis of observed emission lines and derived column 
densities. Due to the complex interplay between chem¬ 
istry, temperature, shielding, and radiative transfer in the 
emitting lines, it cannot be exluded that more massive 
optically thick disk models yield similar column densities 
and line emission, but this is left for future modeling. 

Assuming that CO emission scales linearly over a small 
range of disk mass, the comparison with observations of 
CO and [CI] made for /3 Pictoris constrains the gas mass 
between 0.2 and 4 M®. The question remains: do the re¬ 
sults indeed lead to the conclusion that the j3 Pictoris disk 
is depleted in gas, or is there still some physical input miss¬ 
ing in the models? 

Evidence that our models are still not complete comes 
from the recent observations of Le cavelier des Etangs et 
al. ( |200l[ ) and Deleuil et al. (2001), who have shown that 
P Pictoris might be an active star. The detection of the 
O vi emission doublet at 1035 A and of C n and C ill 
lines with the FUSE satellite, lead to the conclusion that 
P Pictoris may have a chromosphere. Even though the ob¬ 
served O VI lines do not overlap with any of the CO pho¬ 
todissociation lines, the continuum of the chromosphere 
increases the stellar UV flux and hence the amount of CO 
photodissociation. H 2 will be less affected due to its more 
efficient self-shielding. The existence of a chromosphere 
in P Pictoris would increase the photodissociation rate for 
CO as well as the ionization rate for carbon, hence leading 
to lower CO and Cl column densities. 


FUSE is not sensitive enough to detect the continuum 
flux shortward of 1100 A in P Pictoris. This makes it dif¬ 
ficult to estimate the strength of the additional chromo¬ 
spheric UV radiation field. Nevertheless, the upper limit 
for a continuum given by the FUSE spectra (at 1000 A 
~ 3.2 x 10 -11 erg cm -2 s -1 Hz -1 ), could still hide a chro¬ 
mosphere of about 10 times higher than the photospheric 
flux included so far in the models. Section ^ shows that the 
inclusion of an interstellar UV radiation field decreases the 
CO emission already by a factor of 20 (compare Fig. |T^ n 
and c). Hence the presence of a chromosphere in the star 
P Pictoris, which can be according to the FUSE data up 
to a factor 10 in units of the interstellar radiation field, 
may allow larger disk masses to be consistent with the cur¬ 
rent CO J=2-l and 3-2 observational limits. Moreover it 
can enlarge the [Cn] 157.74 /im line significantly and thus 
bring it closer to the observed integrated intensity. 

Heap et al. ( |2000| ) report the detection of an inner warp 
in the disk around P Pictoris. From STIS coronographic 
observations, they deduce a two component disk model: a 
main outer disk and a fainter inner disk, which is inclined 
by 4-5° with respect to the main disk and extends to about 
80 AU. Hence, the inclined inner disk will receive more 
stellar UV radiation, leading to less CO, Cl and more 
Cil. In the case of an inner warp, the UV absorption-line 
studies would miss some of the inner disk material and un¬ 
derestimate the total disk mass. Both effects allow larger 
disk masses to be consistent with the present observations. 

So we are left with the final question: do the dusty 
disks seen around A stars (“Vega-type” stars) still contain 
gas? The previously drawn conclusion that these disks are 
debris disks which contain less gas than dust, is mainly 
based on CO conversion factors appropriate for molecular 
clouds. But the low mass disks around A-type stars re¬ 
semble more a photodissociation region than a molecular 
cloud. From the models presented here and the discussion 
of the various observations for P Pictoris, it is difficult to 
keep the gas-to-dust ratio as high as 100. On the basis of 
the present models, we argue for a gas-to-dust ratio be¬ 
tween 0.5 and 9 for P Pict oris, assuming a dust mass of 
0.44 M® (Chini et al. 1991 ). This ratio depends of course 
on the dust mass and on th e basis of different modeling ap¬ 
proaches Dent et al. ( 2000|) f ound a dust mass of 0.04 M®, 
while Li & Greenberg (1998) obtained 0.33 M®. Any chro¬ 
mosphere included in the models will - depending on 
its strength — raise the gas-to-dust value. Moreover, the 
2 M® model with the interstellar radiation field suggests 
that the inclusion of a chromospheric UV radiation field 
can solve the remaining problem with CO, Cl, and Cil 
simultaneously. This shows the power of using observa¬ 
tions of several species to constrain the gas mass instead 
of using only one, namely CO. 

For Vega, the non-detection of CO gave an upper limit 
for the gas mass of 7 x 10 -3 M® using a CO abundance 
of ~ 10 -4 (Yamaslrita et al. 1993). Our models show, that 
CO is entirely photodissociated in Vega and hence this 
upper limit has to be regarded with caution. Following the 
models presented in this paper, most of the carbon is in 
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form of C+. Hence, we use [Cll] observations to constrain 
the gas mass in the disk around Vega and obtain an upper 
limit of 0.2 M s . This gives and upper limit of 33 for the 
gas-to-dust ratio, assuming a dust mass of 6 x 10 -3 M 0 
(Chini et al. fl990| ) an 
10 -3 M 0 (Dent et al. 

In order to imprcr 
and to constrain the gas mass, more suitable gas tracers 
than CO are needed. While absorption studies only probe 
the disk material in the line of sight, emission line studies 
contain information on the disk structure as a whole and 
are therefore better suited to constrain the disk models. 
The calculations presented in this paper for a 2 M 0 disk 
with a gas-to-dust mass ratio of 100 show that the best 
tracers among those considered here (CO, C, and C + ) are 
C + and C. This result does not depend on the geometry 
of a nearby disk, because the beams of APEX, SOFIA 
and Herschel are small. Moreover the conclusion holds for 
disks around A stars with a moderate radiation field like 
/3 Pictoris as well as for disks that are exposed to a larger 
UV flux like around Vega. 
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